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Measurements of the 1 H-NMR spectrum of a small (~ 4 fj,g) single crystal of the organic conductor 
A-(BETS)2FeCl4 are reported with an applied magnetic field Bo = 9 T parallel to the a-axis in the 
ac-plane over a temperature (T) range 2.0 — 180 K. They provide the distribution of the static 
local magnetic field at the proton sites in the paramagnetic metal (PM) and antiferromagnetic 
insulator (AFI) phases, along with the changes that occur at the PM— AFI phase transition. The 
spectra have six main peaks that are significantly broadened and shifted at low T. The origin of 
these features is attributed to the large dipolar field from the 3d Fe 3+ ion moments (spin Sd = 
5/2). Their amplitude and T— dependence are modeled using a modified Brillouin function that 
includes a mean field approximation for the total exchange interaction (Jo) between one Fe 3+ ion 
and its two nearest neighbors. A good fit is obtained using Jo = — 1.7 K. At temperatures below 
the PM— AFI transition temperature Tmi = 3.5 K, an extra peak appears on the high frequency 
side of the spectrum and the details of the spectrum become smeared. Also, the rms linewidth 
and the frequency shift of the spectral distribution are discontinuous, consistent with the transition 
being first-order. These measurements verify that the dominant local magnetic field contribution is 
from the Fe 3+ ions and indicate that there is a significant change in the static local magnetic field 
distribution at the proton sites on traversing the PM to AFI phase transition. 

PACS numbers: 75.30.Kz, 75.50.Ee, 76.60.-k, 71.30.+h 



I. INTRODUCTION 

The organic conductor A-(BETS) 2 FeCl 4 , where BETS 
is bis(ethylenedithio)tetraselenafulvalene (CioS4Se4H 8 ), 
is of considerable interest because of the properties 
related to the coexistence of the large magnetic 3d 
Fe 3+ moments (spin 5d=5/2) of the inorganic anions 
(FeCl 4- ) and the conduction 7r-electrons (spin 5^=1/2) 
in the donor molecules from the BETS. 1 ' 2 ' 3 ' 4 ' 5 ' 6 It has 
an unusual phase diagram, including an antiferromag- 
netic insulating (AFI) phase, a paramagnetic metallic 
(PM) phase, and a field-induced superconducting (FISC) 
phasei^ Also, results interpreted in terms of a ferroelec- 
tric phase transition in the metallic phased at 70 K and a 
relaxor ferrolectric behavior— at 30 K have been reported. 

A mechanism used to explain the FISC phase below 
5 K in A-(BETS)2FeCl4 is based upon the Jaccarino- 
Peter (J-P) compensation effect 10 operating in a two- 
dimensional (2D) systemiSiiiiiS In this model, the neg- 
ative exchange interaction (J^d) between the paramag- 
netic 3d Fe 3+ moment (g/isSd) (g is the Lande factor and 
\xb is the Bohr magneton) and the conduction 7r-electrons 
in the BETS molecule (it— d interaction) 1 ^ generates a 
large magnetic field (B^d) that cancels most of the exter- 
nally applied magnetic field (Bo) when the latter is large 
(Bo ^17 — 45 T) and aligned parallel to the ac plane. 
This suggests that B^d (-B^d = J-nd < Sd >/g^B, where 
< Sd > is the average value of the Fe 3+ spin polariza- 
tion) at the conducting 7r-electrons is on the order of ~ 
30 T and its direction is antiparallel to Bo- However, this 
interpretation has been challenged because the supercon- 



ducting state in A-(BETS)2FeCLi can be destroyed by a 
very small amount of out-of-plane magnetic field (~ 0.1 
T)i£ and the alignment of the paramagnetic Fe 3+ mo- 
ment follows closely that of Bo- It has also been proposed 
that the Larkin-Ovchinkov-Fulde-Farrell (LOFF) phase 
is present near the boundary of the FISC phase. 13 - 14 

The magnetic PM— AFI phase transition, which occurs 
at Bq < 11 T (Neel temperature Tjv), coincides with a 
metal-insulator (MI) transitionii 5 * 1 ^ (transition temper- 
ature Tmi)- The property Tjv = Tmi indicates that the 
MI and AFI transitions are cooperative transitions! 1 * 1 ^ 6 - 
Thus, it is expected this PM— AFI transition is also 
a result of the tt— d interaction, since the study of its 
iso-structural nonmagnetic and non-3d-electron analog 
A-(BETS) 2 GaCl 4 shows that it exhibits a behavioriLiS 
that is completely different from that of A-(BETS)2FeCl4. 
However, the detailed role of the 7r— d interaction for the 
PM— AFI phase transition is not yet clear because it has 
not yet been established by direct experimental evidence. 

Prior reports of proton NMR measurements on A- 
(BETS^FeCLt include the spectra of a large (~ 6.5 mg) 
aggregate of crystals aligned along the c-axis in a mag- 
netic field 7 of 2.2 T and a preliminary report^ of our 
measurements on a ~ 4 /xg single crystal, for which the 
spectrum results are presented and analyzed in detail 
here. A comparison of the results from the aggregate 
sample 7 and the significantly different ones on the single 
crystal reported here is discussed is Section IV. F. 

The crystal structure^ 7 , of A-(BETS) 2 FeCl 4 , shown 
in Fig. 1, is triclinic with the space group PI and the 
lattice constants at 298 K a = 16.164(3) A, b = 18.538(3) 
A, c = 6.5928(8) A, a = 98.40(1)°, (3 = 96.67(1)°, 7 




FIG. 1: (color online) (a) The crystal structure of A- 
(BETS^FeCU with all 32 hydrogen atoms labeled in a unit 
cell. Atoms in color: Red — Fe; Green —CI; Blue — H; Grey 
-C; Yellow -S; Pink -Se. (b) The BETS molecule, which 
identifies the Se, C, and S atoms. 

= 112.52(1)° and V = 1773.0(5) A 3 . There are four 
BETS molecules and two Fe 3+ ions per unit cell and the 
BETS molelcules are stacked along the a- and c-axes to 
form a quasi-stacking fourfold structure. The conducting 
layers, comprised of BETS, are sandwiched along the b- 
axis by the insulating layers of FeCl 4- anions. The least 
conducting axis is 6, ac is the conducting plane, and the 
easy axis of the antiferromagnetic spin structure is ~ 30° 
away from the c axis (|| needle axis of the crystal), 20 

In this paper, measurements of the proton NMR spec- 
trum of a single ~ 4 ^g crystal of A-(BETS)2FeCl4 are 
reported with a magnetic field Bo = 9 T applied paral- 
lel to the a-axis in the ac-plane over a temperature (T) 
range 2.0—180 K. Along with our preliminary report, it 
is the first published report of proton NMR in a single 
crystal of this material. It is a challanging experiment 
involving an extremely small sample size (~ 0.018 mm 
x 0.065 mm x 1.2 mm). These measurements probe the 
distribution and the origin of the static local magnetic 
field at the proton sites in the PM and AFI states as well 
as across the PM— AFI phase transition. The observed 
properties should help to establish a microscopic model 
for the PM— AFI phase transition. 

One important result of this investigation is that the 
dominant local magnetic field at the proton sites comes 
from the large dipolar field of the 3d Fe 3+ ion moments. 
A mean field model based on the dipolar field of the Fe 3+ 
moments is presented and used to calculate the proton 
NMR spectrum. It provides a good fit to the measured 



spectra. Besides this, the total exchange constant Jo be- 
tween an Fe 3+ ion and its two nearest neighbors is deter- 
mined to be Jo ~ — 1-7 K from a fit to the spectrum data. 
These measurements also show that there is a significant 
change in the static local magnetic field distribution at 
the proton sites across the PM— AFI phase transition. No 
proton NMR evidence of a ferroelectric phase transition 
at 70 K is observed in these measurements. 

The rest of this paper is organized as follows. Section 
II presents the experimental details and Section III has 
the experimental results for the proton NMR spectra, 
including frequency distributions, shifts, and linewidths. 
Section IV presents the model for the spectrum, along 
with the comparison with the measured spectra. The 
conclusions are stated in Section V. 



II. EXPERIMENTAL DETAILS 

The needle- like single crystal A-(BETS)2FeCl4 samples 
were prepared as described by Montgomery et al. with 
a standard electrochemical oxidation method^ The di- 
mensions of the sample used for these 1 H-NMR measure- 
ments are (1.2 ± 0.1) mm x (0.065 ± 0.010) mm x (0.018 
± 0.005) mm, which corresponds to (3.8 ± 1.8) fig in mass 
and (2.7 ± 1.3) x 10 16 protons. 

The NMR coil used was 40 turns of 0.025 mm diame- 
ter bare copper wire wound on a 0.075 mm diameter wire 
form. The coil was held to the rest of the probe circuit 
by two 125 fim diameter Cu wire leads with Teflon in- 
sulation, which was removed close to the ends where the 
coil was soldered to them. Commercial pure acetone was 
used for cleaning the coil and its surroundings when the 
NMR coil was set on the probe to reduce the spurious 
proton signals relative to the signals from such a small 
sample. 

Finally, a single piece of the needle shape A- 
(BETS^FeCLj single crystal sample was slid into the coil, 
carefully aligned close to Bq || a in the ac-plane, and held 
in place with a very small amount of commercial Apiezon 
grease on each end. The orientation is done visually un- 
der a microscope with an estimated uncertainty of ~ ± 
5°. [Note: the ac-plane is ~ in the sample surface plane 
which has the largest surface area, the c-axis 1 1 the needle 
direction of the sample, and the angle between the a- and 
c-axes is (3 (/3 = 96.67(1°) at 298 K)]£iL2I 

As shown in the preliminary report^ the spurious pro- 
ton signal was estimated to be less than ~ 4% of the sig- 
nal from the A-(BETS)2FeCi4 sample by comparing the 
signal with and without the sample in the coil. Thus, 
in these measurements the spurious proton signal has an 
insignificant size. 

The 1 H-NMR frequency-swept spectra were obtained 
using standard spin-echo techniques carried out with a 
spectrometer and probe built at UCLA. Since the pro- 
ton has a gyromagnetic ratio ji = 42.5759 MHz/T, the 
frequency v for the excitation pulses used for the spec- 
trometer is near v <~ vq = 7/ Bo = 382.6935 (MHz), where 
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is the proton Larmor frequency in the external field. 
The value of Bq used in this experiment was Bo — 8.9885 
T (for simplicity, often referred to here as 9 T). 

Because the NMR spectrum covers a wide range in fre- 
quency up to 14 MHz (3.3 kG), short rf pulses and a wide 
receiver bandwidth (± 1 MHz) were used to record the 
spin-echo signals. The pulse sequence that optimized the 
height of the spin echo used to record the NMR signal 
was a 0.2 /is tt/2 pulse (B 1 = 294 G, 1.25 MHz proton 
frequency) followed by a 0.3 /xs pulse separated by a time 
interval r (r ~ 5 /is) for most of the measurements. For 
a viable signal-to-noise ratio, each echo signal was aver- 
aged 2000 times at 180 K and 128 times at 4.2 K and 
lower temperatures. The uncertainty associated with the 
signal-to-noise ratio is probably the main source of error 
in the data. The uncertainties include ~ ± 1% in T and 
± 5° in the field alignment. 

At low T, the spectrum is very wide (~ 12 MHz) and 
the frequency sweep covered a range as high as 370 to 400 
MHz and used a typical frequency step for each acqui- 
sition of 0.2 — 0.5 MHz. When a wide frequency sweep 
range was used, the probe circuit was retuned every 4 
MHz to maintain a uniformly high sensitivity (above 
85%) for recording the proton spectrum. The spec- 
tra were analyzed with frequency-shifted and -summed 
Fourier transform processing^ 

III. RESULTS 

A. 'H-NMR spectra 

Figure 2 shows the normalized 1 H-NMR absorption 
spectra [x (f)] of a A-(BETS)2FeCi4 single crystal as a 
function of the frequency shift v — Vq with Bo = 8.9885 
T parallel to the a-axis in the ac-plane for T = 180 K to 
2 K. 

Over most of this range, the spectra have six main 
peaks which can be divided into two groups (low fre- 
quency side and high frequency side) with 3 main peaks 
for each. As T is lowered from 180 to 4.2 K in the 
PM state, the spectrum broadens significantly, its cen- 
ter shifts to lower frequency and the splitting between 
the peaks increases. For such a complex spectrum, a rea- 
sonable measure of its center is the first moment, < v >, 
given by^ 3 - 

Ei {Vi) 

where i indexes equally spaced frequency steps. The av- 
erage shift of the spectrum (Ai/) indicated by the dashed 
line in Fig. 2, is Au =< v > — v§. Also, some additional 
weak structures gradually develop at lower T. As the 
sample is cooled further into the AFI phase (T < Tmi 
= 3.5 K), the details of the spectrum become somewhat 
smeared, an additional peak at 4.25 MHz appears and it 
grows larger with further cooling (solid line in Fig. 2). 
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FIG. 2: Normalized H-NMR absorption spectra of a single 
crystal of A-(BETS)2FeCLt as a function of v — va (vo = 
382.6935 MHz) from 180 K to 2 K with B = 8.9885 T parallel 
to the a-axis in the ac-plane. The solid vertical line at v — v o 
= 4.25 MHz indicates a new peak below T < Tmi = 3.5 K, 
and the dashed line shows < v > for each spectrum. 



As discussed in later sections, this relatively complex 
spectrum is caused mainly by the dipolar field of the 
3d Fe 3+ ion electron spin moments (S^ = 5/2, g m 2) 
at 16 magnetically inequivalent proton sites in both the 
PM and AFI phases. In both phases, these moments 
are present, the major differences being that they should 
have long-range order and a different orientation in the 
AFI state. 

At the lowest T in the PM state, the average mganeti- 
zation of each Fe 3+ ion is almost completely saturated by 
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Bq. Therefore, the main changes in the NMR spectrum 
associated with the AFI phase are caused by changes in 
the orientation of the Fe 3+ average spin moments, not by 
changes in their magnitude. 

In the AFI phase (T < Tmi), the changes of the spec- 
trum vs T, such as the smearing of the spectrum details 
and the growth of the new peak at the high frequency 
side, occur gradually rather than suddenly. This behav- 
ior indicates a continuous change of the local magnetic 
field distribution at the proton sites as T is lowered. 
Its origin is probably the development of the canted- 
antifcrromagnetic phase as observed from the field depen- 
dence of magnetization and capacitance measurements 
reported 2 ^ for A-(BETS) 2 FeCl 4 . 



B. 1 H-NMR frequency shift 

Figure 3 shows the frequency shift magnitude \Av\ as 
a function of T, and the inset of Fig. 3 plots l/\Av\ vs 
T. A clear discontinuity of Av at T = Tmi is seen. The 
error bars are our best estimate of the uncertainty in our 
data analysis. 

The NMR frequency at the «th proton site is 77 x |Bj| , 
where is the total magnetic field at the site. Since B; 
is the sum of the large applied Bo plus the much smaller 
field generated by the sample (AB ), it is easily shown 
that the field shift caused by the magnetic properties of 
the sample is the component of its field parallel to Bo 
(A£>||). The values of Av in Fig. 3 represent the average 
of A-Bm over the proton sites (< AB11 >) with 



N 

X 



■a 
a 



1.0 



S3 



>, 
u 
a 

<u 

a 
cr 

<u 
u 



z 

X 



0.1 



X - (BETS) 2 FeCl 4 



B =9T, II a 



• experimental lAvl 
— fit with J =-1.7K 



1.2 



S 1.0 



0.8 



U 0.6 



0.4 




Ji-(BETS) 2 FeCl 
B =9T, II a 



T = 3.5 K 

MI 

I ..«»-""* 
J 



5 10 IS 20 
temperature T (K) 



f 



Av = 7/ < A£?|| > 



(2) 



1 10 100 200 

temperature T (K) 

FIG. 3: (color online) T— dependence of \Au\ for a single crys- 
tal of A-(BETS) 2 FeCl 4 from 180 K to 2 K with B = 8.9885 
T parallel to the a-axis in the ac-plane. The solid (red) line 
is a fit to Eqs. (6) — (10) with Jo — — 1.7 K. The inset shows 
the discontinuity of 1/|Ai4 vs T across the PM— AFI phase 
transition at 3.5 K. The dashed (blue) line is a guide to the 
eye. 



Thus, the origin of Av is the distribution of A_Bii over 
the proton sites. 

At 180 K, Av has a value of -(0.13 ± 0.02) MHz, while 
at 4 K it reaches -(1.750 ± 0.005) MHz. Below T M i = 
3.5 K, Av has a rather weak temperature dependence, 
as shown in the inset of Fig. 3. The sudden decrease 
of 1/|A^| at Tmi indicates an increase of the average 
static local magnetic field at the proton sites due to the 
PM— AFI phase transition. Its sharpness is evidence that 
the PM— AFI phase transition is first order. 

The negative sign of Av indicates that the direction of 
the average static local magnetic field is opposite to the 
external magnetic field. But this does not mean that the 
average static moment is negative. The negative sign is 
determined by the proton positions relative to the Fe 3+ 
ions in the crystal lattice. 



C. 'H-NMR spectrum linewidth 

A reasonable quantitative measure of the width of the 
local field distribution that generates the proton spec- 
trum of A-(BETS)2FeCi4 is the root mean square (rms) 



linewidth, A/ rms , i.e. the square root of the second mo- 
ment, < Av 2 > : / 2 , given by 23 ' 24 



A/ rms ^< A^ 2 >^- 



E>j- < " >) 2 X>j) 

YriX"(Vi) 



1/2 



(3) 



The measurements of A/ rms for protons as a function 
of T in A-(BETS) 2 FeCl4 are shown in Fig. 4, where 
l/A/ rms is plotted as a function of T. At 180 K, A/ rms 
has the value of (0.205 ± 0.005) MHz and it reaches (2.50 
± 0.10) MHz at T M i = 3.5 K. The error bar is estimated 
to be equal to or smaller than the size of the marker 
below - 50 K. 

Across Tmi = 3.5 K, as for \Av\ in Fig. 3, there is a 
sudden increase of A/ rms in Fig. 4 that also indicates a 
first order phase transition. 

There are several important properties of the data 
shown in Figs. 2—4. At T > 50 K, it can be shown 
that the T— dependence of all of them (< v >, \Av\, 
and A/ rms ) follows the Curie- Weiss relation with a Curie- 
Weiss temperature O ~ 5.5 K. They also follow the Bril- 
louin function [_Bj(x)] behavior using the parameters for 
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Fe 3+ , and Bj(x Q ) is 
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FIG. 4: (color online) T— dependence of l/A/ rms for a single 
crystal of A-(BETS) 2 FeCl 4 from 180 K to 2 K with B = 
8.9885 T parallel to the a-axis in the ac-plane. The solid 
(black) line shows the best fit to the data from Eqs. (6) — (10) 
with Jo = — 1.7 K. The inset shows the low T data across 
the PM— AFI phase transition, with a dashed line as a guide 
to the eye. 
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The effect of Jo for an antiferromagnetic exchange inter- 
action between the nearest neighbor Fe 3+ moments can 
be modeled as an additional magnetic field component 
(£?') antiparallcl to B r& given by 



o/ \Jo\JkB d / \ \Jo\Jk B 

B = Bj(x) ss 

3Mb gm 



Bj{x ), 



(7) 



where 



x 



Jgfi B (Bo-B') 
k B T 



(8) 



and the right side of Eq. 7 has been used for B'. Because 
the value obtained later for B' is substantially smaller 
than J? , this approximation is a reasonable one. Thus, 
the mean field modified Brillouin function used to model 
our data is Bjm (x) and the corresponding formula to fit 



A/rms = C" M(x) = C'N Agf M B J B JM (x), 



(9) 



where 



the Fe 3+ spins. This result is a strong indication that the 
dominant contribution to ABo is the dipolar field at the 
proton sites from the magnetization of the Fe 3+ spins. 

An important issue for characterizing the data is that 
they extend well into the low T regime; i.e., at 9 T, the 
Zeeman splitting between the highest and lowest Fe 3+ 
spin states is 5 x g[i B B^jk B — 60.5 K, where k B is the 
Boltzmann constant. For this reason, instead of the more 
commonly used Curie— Weiss law, which is a high T ap- 
proximation (T >> gfJ- B Sd/k B ~ 30.2 K), we character- 
ize the data with the Brillouin function, which includes 
magnetic saturation at low T. Also, since there is a sig- 
nificant total antiferromagnetic (negative) exchange in- 
teraction Jo among the Fe 3+ moments^^ the Brillouin 
function used here is modified to include a simple, ap- 
proximate mean field correction. 

For non-interacting moments, the magnetization 
[M(x )] is given by2& 

M(x Q ) = N A9 fi B J Bj(x ), (4) 

where Na is Avogadro's number, J = Sd = 5/2 for the 



Bjm(x) =Bj(x +x'), (10) 

and C and Jo are adjusted to give the best fit to the 
data. 

As shown in Fig. 4, from the T— dependence of A/ rms 
the best fit is obtained with J = -(1.7 ± 0.2) K. This 
corresponds to a maxmium total exchange field of ~ — 3 
T below ~ 5 K. The negative sign of Jo indicates that the 
Fe 3+ ions have an antiferromagnetic (AF) exchange inter- 
actions with their nearest neighbors. Since each Fe 3+ ion 
has two nearest neighbors^ (actually three closest ones: 
two are at ~ 6.6 A away, and one at ~ 7.6 A away) the 
exchange constant Jdd between each pair of Fe 3+ ions 
is, Jdd ~ ^o/2 = —(0.85 ± 0.10) K, which agrees with 
theoretical expectations (~ — 0.64 K). 2 - The parame- 
ter C has a fitted value of (100 ± 6) [(mol.Fe/emu)Hz] 
[Note: The units for M(x) are emu/mol.Fe; see Eqs. (4) 
and (9)]. The overall difference between the fit and the 
A/ rms data is below ~ 5%, except near the phase tran- 
sition, where it is ~ 10%. 

Similarly, \Ai>\ (Fig. 3) is well characterized by Eqs. 
(6) -(10) from 180 K down to 10 K with the same fit 
value of Jo- But its deviation is slighly larger below ~ 10 
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K. This is possibly caused by not including the demag- 
nitization and Lorentz fields to the local field B' in Eq. 
(7) for the shift. 

The property A/ rms > |Af| shows that there is 
a broad static local magnetic field distribution in A- 
(BETS^FeCLj. It occurs because there are 16 inequiv- 
alent 1 H-sites at which the dipolar field from the Fe 3+ 
moments has a large variation. 

As discussed in more detail below, these proton shift 
properties support the conclusion that they are domi- 
nated by the dipolar field from the 3d Fe 3+ ion electron 
moments. The sudden change in the spectrum at Tmi 
= 3.5 K, as well as seen from those proton shift proper- 
ties, reflects a comprehensive change of the static local 
magnetic field distribution at the 1 H-sites due to the AF 
ordering of the Fe 3+ electron spins. 



IV. DISCUSSION 

In this section, the local magnetic fields at the proton 
sites and the 1 H-NMR spectrum are calculated, and the 
nature of the AFI phase transition is discussed. 



A. Model for the 1 H-NMR spectra 

Generally NMR spectra are determined by the values 
of static local magnetic field and the distribution of the 
local magnetic field at the nucleus sites in the studied ma- 
terial. Here, a model for the spectrum is presented and 
applied to A-(BETS)2FeCLi. It considers all possible ma- 
jor sources which include the dipolar field of the Fe 3+ , the 
exchange interactions with the Fe 3+ ion and 7r-electrons, 
the dipolar field of the neighboring proton nuclei, and 
the demagnetization and Lorentz contributions^SiS&iSl 
to the local field at the proton sites. 

In A-(BETS) 2 FeCl 4 single crystals^ there are 16 in- 
equivalent proton sites per unit cell (see Fig. 1). Thus, 
up to 16 different lines in 1 H-NMR spectrum can be ex- 
pected. Each of these 16 protons will have, in general, 
a different shift in the NMR frequency depending on its 
position in the crystal lattice. 

The Hamiltonian Hi of the system for the 1 H-NMR 
can be expressed as^ 3 - 



Hi 



Hiz+Hii+Hf*- 



H Id- 



H I1T +H 



dcm 



H 



Lor 



(ii) 



where Hiz is the Zeeman Hamiltonian of the *H nuclei in 
Bo , Hn is the proton-proton nuclear dipolar interaction 
Hamiltonian, Hj^ and ijj^ are the dipolar coupling and 
transferred hyperfine coupling from the 3d Fe 3+ electrons 
to the protons, respectively, Hi^ is the hyperfine coupling 
of the proton nucleus to the BETS 7r-electrons, and the 
last two terms, H dcm and H Lor , are the bulk demagneti- 
zation and Lorentz contributions, respectively^ 3 *^ All of 
these terms contribute to the static local magnetic field 
at the proton sites and all but the first cause the 1 H-NMR 
frequency shifts. 



Because of the small atomic number of the hydrogen 
nucleus (Z — l)p2, it is expected that the the proton 
hyperfine couplings to the 7r-electrons (Hi^) and to the 
Fe 3+ electrons (-ff/j) are negligible. For this reason, Hi 
is simplified to 



Hi » Hiz + Hu + H 



dip 
Id 



H 



dcm 



H 



Lor 



(12) 



Among these terms it is expected that the dipolar field 
of the Fe 3+ ion electron spins, i.e. the contribution of 
Hamiltonian Hf 1 ^ , to be the dominant source contribut- 
ing to the static local magnetic field at the proton sites, 
as seen from the analysis in the following sections. 

Note that any interactions, both direct and indirect, 
between the 7r-electrons and the Fe 3+ ions (ir—d interac- 
tion) or between the Fe 3+ ions (d— d interaction) will 
affect the polarization of the Fe 3+ electron moments, 
therebye modifying the dipolar field from the 3d Fe 3+ 
ions at the proton sites. In what follows, the effects of 
these interactions are considered, and the total d— d ex- 
change interactions including those2& through the Cl~ 
and the conduction it— electrons (RKKY interaction) are 
included in the calculation using a mean field approxi- 
mation to modify the Brillouin function. Since the di- 
rect 7T— d interaction is considered to be small because 
the magnetization of the conduction tt— electrons is very 
small compared to that of the Fe 3+ , it will not be in- 
cluded. 

The dipole moment flj of the Fe 3+ ion j produces a 
magnetic field Bij at the proton site i given by2& 



Bij - 



(13) 



where is the position vector from the proton site i to 
the Fe 3+ ion site j. 

Thus, the total dipolar field < Bi > at the proton site 



< Bi> = ^< B^ >, 



E< 

3 



3fy(/I-ry) fa 



(14) 



>, (15) 



and the magnetization M(x) of the Fe 3+ moments is 

£ 3 - < N > 



M(x) 



V 



(16) 



By considering Eqs. (6) — (10) and (13) — (16), and assum- 
ing that the magnetization M(x) has the same direction 
as Bo, one obtains 



< Bi > r 



gfi B J Bjm(x) 

Bn 



+N 

£ 

3 = -N 



3rij(B ■ fjj) _ do 

3 



Ta 



(17) 
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where i = 1, 2, 16, which indexes the proton sites, and 
a; is a T— dependent variable as defined by Eqs. (6) — (8). 

Equation (17) is used to calculate the local field at all 
T in the PM phase. The mean field approximation for 
the Fe 3+ exchange interactions is included in Bjm(x). 

At high T (gn B B J « k B T), Eq. (17) can also be 
written as 



< Bi >? 



(g^) 2 J(J+1) 

3fc s r+e 



+N 

£ 

j=-N 



3fij(B -rij) B t 



(18) 

where is the Curie- Weiss temperature. 

The field Bo can be expressed as Bo = £>o(sin 8 cos (f>i + 
sin 8 sin </>j+cos 8k), where 8 and tj> are standard spherical 
coordinates in the Cartesian system^ and the dipolar 
field components < Bi > x , < Bi > y , and < Bi > z along 
x, y, z directions, respectively can be calculated for each 
of the 16 inequivalent protons sites with Eq. (17) by 
considering the coordinates of all the inequivalent proton 
sites and the Fe 3+ ion positions. 

Since the values of < Bi > obtained in the next section 
obey | < Bi > \ « B , the contribution of < Bi > to 
the shift of the proton spectrum comes only from the 
component of < Bi > || B (Bf, lp ). 



B. Calculated local magnetic field at the proton 
sites 

Figure 5 shows the calculated component of the Fe 3+ 
ion dipolar field parallel to Bo at the 16 inequivalent 
proton sites HI, H2, H16 in the crystal lattice using 
Eq. (17) at T = 20 K and B = 8.9885 T. The number 
of unit cells included is (2x100 + l) 3 ; i.e., N = 100. The 
crystal ac— plane is chosen to be in the xz— plane in the 
transformation of the lattice triclinic coordinates to the 
Cartesian coordinates used for the calculation. 

The values of the components of the Lorentz field, 
B^ m , and the demagnetization field, B|j oln , parallel to 
B , based on the shape of the needle-shape single crys- 
tal sample and the magnetization M(x) [Eqs. (5) — (10)], 
are +108.3 G and —32.5 G, respectively 28 ' 31 as shown 
in Fig. 5. The net shift from both of these contribu- 
tions is their sum, i.e. +75.8 G, which is small, but not 
completely negligible. The small spacial variations^ of 
the demagnetization field across the sample have been 
neglected. 

Also, the calculated contribution of Hu is only < 3 G 
among the 16 inequivalent proton sites, as confirmed by 
the spin-echo decay measurements Since it is so small, 
it is also neglected here. 

Thus, B^ lp obtained from Eq. (17) (see Fig. 5) is the 

dominant contribution to the structure of x"iy)- Since 
B^° T and B^ cm are nearly constant over all the proton 
sites, they have a negligible effect on the structure of 
x" their contribution constitutes a shift in frequency 
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FIG. 5: (color online) Calculated dipolar (Bjj lp ), Lorentz 
(B^° T ), and demagnetization field (_B|j cm ) components par- 
allel to Bo at each of the 16 inequivalent proton sites (index 
i) in a single crystal of A-(BETS) 2 FeCl 4 at 20 K. The align- 
ment of Bo = 8.9885 T is parallel to the a-axis in the ac-plane. 
The solid and dashed lines are guides to the eye. 



or field but with almost the same amount for each proton 
sites (~ 25% to the average Av at 20 K). 

The values of -B|j lp at 20 K cover a wide range of field, 
from ± 20 G up to ± 600 G, depending on the proton po- 
sitions in the crystal lattice. For example, at the proton 
site H15 it is ~ —600 G, and ~ +500 G at the proton site 
H12. The positions of these proton sites are all labeled in 
Fig. 1. The proton site (H15) that has the largest B^ lp 

from the Fe 3+ ions is that which is closest to the nearby 
Fe 3+ ion plane in the crystal lattice. The range of field 

(~ —600 h600 G) that B|j lp covers corresponds to a 

range of frequency of ~ 5 MHz. 



Calculated 1 H-NMR spectra from the dipolar 
field contributions 



In this section the model for the proton absorption 
spectrum [XmodM] i s calculated. The first step is to 
calculate the static local magnetic field B^ 1P at each of 
the 16 inequivalent proton sites, as described in Sections 
IV A and B. The second step is to convolve this field 
distribution with a set of Gaussian functions, yi, each 
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FIG. 7: (color online) Comparison of the calculated Xmod(^) 
[upper (blue) smooth line] with the measured 1 H-NMR ab- 
sorption spectrum \ i v ) [lower (red) line] in a single crystal 

of A-(BETS) 2 FeCi4 at 20 K. For XmodO), tne val ue of B 
is 8.9885 T aligned with 9 = +85° and <j> = +5°, which is 
close to the a-axis and near the ac-plane. The vertical arrows 
indicate the center (first moment < v >) of each line. 



These angles are selected because they are close to what 
is needed for comparison with the measurements. As can 
be seen in Fig. 6, X m od( u ) ^ s fairly sensitive to 9 and <fi 
which are determined by the direction of Bq . 



FIG. 6: (color online) Calculated proton XmodM from the 
3d Fe 3+ ion electron dipolar contributions in single crystal 
A-(BETS) 2 FeCl 4 at 20 K with B = 8.9885 T aligned close 
to the a-axis (9 = 90° ±5°) and near the ac-plane ( — 15° < 
4> < +20°). The Larmor frequency v = 382.6935 MHz. 



of which has a maximum amplitude of 1 and the same 
width 8 at each proton site In this case, 



^exp 



i=l 



2S 2 



(19) 



Plots of x mo dM! calculated with this model and the 
value of 5 — 0.15 MHz, are shown in Fig. 6 for T = 20 
K with Bq — 8.9885 T aligned close to the a-axis (9 = 
90° ± 5°) and near the ac-plane (-15° < (j) < +20°). 



D. Comparison of the model to the measured 
spectra 

Figure 7 shows a comparison of the calculated X mo d(^) 
(upper) using B = 8.9885 T with 9 = +85° and <f> = +5°, 
and the measured x i v ) (lower) at 20 K for = 382.635 
MHz. The calculated result includes the dipolar field 
from the Fe 3+ electrons, and corrections for the Lorentz 
and bulk demagnetization fields. 

There are several small differences between the mea- 
sured and calculated model spectra shown in Fig. 7. One 
is that at 20 K, the shift in the center of the spectrum dif- 
fers by ~ 42 kHz (~ 100 G). This is quite small compared 
to the separation of the outer peaks, which is 4.65 MHz 
(1.09 kG) for both the measured and calculated spectra. 
On the other hand, A/ rms = 1.126 MHz and 1.387 MHz 
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for the model and measured values, respectively, corre- 
sponding to a difference of ~ 18%. Another small dif- 
ference is in the shape of the spectra. These similarities 
in the model and measured spectra strongly support the 
conclusion that the proton spectrum in A-(BETS)2FeCl4 
is dominated by the dipole field of the Fe 3+ ions. 

An important question is what is responsible for the 
small difference between the spectra calculated with the 
model that has been used and the measured spectra. It 
is probably caused by not including all the details of the 
exchange interactions (including the d— d exchange inter- 
actions through the Cl~ and through the BETS conduc- 
tion 7T— electrons) between the Fe 3+ electron spins when 
using the mean field M(x) or B JM (x) to model the Fe 3+ 
electron spin polarization. 

Also, an important aspect of the 7r— d interaction which 
is also part of the d— d exchange interactions is that it 
should be responsible for the coordination of the occur- 
rence of the PM— AFI phase transition. 

It is expected that if all of the exchange interactions 
are included in detail, they will cause a small change in 
the moment of the Fe 3+ ions and their polarization will 
not be precisely along Bo, even in the PM phase. It 
is likely that these additional factors are responsible for 
the relatively small difference between the measured and 
calculated shift and width in the 1 H-NMR spectrum. To 
include them in a model is beyond the scope of this paper. 



E. The local magnetic field in the AFI state and 
the nature of the PM AFI phase transition 

As in the PM state, the local magnetic field at the 
proton sites in the AFI state (T < Tmi — 3.5 K) is also 
dominated by the dipolar field from the Fe 3+ ions, even 
though Eq. (17) is not specific to and may not apply in 
the AFI state. In both phases, the Fe 3+ ion electron spin 
moments are present and the general Hamiltonian Hi of 
the 1 H-NMR system Eq. (11) or (12) applies in both the 
PM and AFI states. 

However, the change of the 1 H-NMR spectra in A- 
(BETS) 2 FeCl 4 is significant. On cooling from the PM 
to the AFI state, the spectrum broadens, the splittings 
are smeared out, and a new peak appears on the high 
frequency side. 

It is unlikely that these changes are caused by varia- 
tions in the Lorentz or demagnetization field because the 
contribution from each of them is essentially the same at 
each proton site, even though they are proportional to 
the Fe 3+ ion electron susceptibility. 

Therefore the changes of the details in the proton spec- 
tra at any T (including in the PM and AFI phases) come 
mainly from the change of the dipolar field of the 3d 
Fe 3+ ions, which in part are changed by the effect of 
7T— d and d— d interactions on the polarization of the Fe 3+ 
moments. The major difference for the Fe 3+ moments is 
that they should have long-range order in the AFI state, 
which is formed by the tt— d and d— d interactionsii^*i& 



The evidence from the change of the 1 H-NMR spec- 
tra, as well as that reflected by the discontinuities of the 
frequency shift Av and the rms linewidth A/ rms , is in- 
dicative of a first order nature for the PM— AFI phase 
transition in A-(BETS) 2 FeCl 4 ^ 3 . 



F. Comparison with other results 

In this section, we first compare our results on a single 
~ 4 /ig crystal with those reported by Endo et ali on a 
large (~ 6.5 mg) aggregate of crystals aligned along the c- 
axis. One major difference is that since their Bo = 2.2 T 
was applied perpendicular to the c-axis, the correspond- 
ing 1 H-NMR spectrum is the average over all alignments 
of Bo perpendicular to the c-axis. It is expected that 
in comparison to the measurements reported here on a 
single crystal, this average will smear out some of the de- 
tails of the spectrum and will generate a broader range 
of Lorentz and demagnetization fields that also have an 
effect on the spectrum. Another important difference is 
that because Bo in our measurements is 4.1 times larger 
than that used in their work, the splittings and shifts of 
the spectral features is expected to be 4.1 times larger; 
i.e., the spectral resolution of our 9 T measurements is 
significantly higher. 

Their results^ showed the onset of a splitting of 
the spectrum near 70 K that was considered an 
anomaly and interpreted as possible evidence for a 
charge disproportionation of the 7r-electrons associated 
with a ferroelectric-type phase transition. They further 
suggested^ that this splitting is from the hyperfine field 
between the proton nucleus and tt electrons and not the 
dipole field of the Fe 3+ ions. 

Our results disagree with their results and interpreta- 
tion in several ways. First, our measurements (Fig. 2) 
show a continuous increase in the splitting and the shift 
(Fig. 3) of the spectrum as T is decreased from 180 to ~ 
5 K. No evidence for an anomaly near 70 K is observed 
and the main features of the spectra are well explained 
in terms of the model based upon the dipole field of the 
Fe 3+ ions. 

Another problem with the interpretation by Endo et 
ali is that if the spectral splitting is caused by the hy- 
perfine field between the proton nucleus and tt electrons, 
it would require a large hyperfine shift of ~ 50 G, corre- 
sponding to a Knight shift of ~ 0.23%. It is very unlikely 
that such a large shift could occur for the very light nu- 
cleus of a hydrogen atom. 

Experimetal evidence that such a large hyperfine field 
does not occur at the proton sites is also indicated by 
comparing the proton spin-lattice relaxation rate (1/Ti) 
at 100 K inia A-(BETS) 2 FeCl 4 and 3 - 4 A-(BETS) 2 GaCl 4 . 
Both materials have the same structure and related prop- 
erties, but without the magnetization of the Fe 3+ ions 
in the Ga compound. If the main fluctuating magnetic 
field at the proton site were from the tt electrons, one 
would expect approximately the same value of 1/Ti for 
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both materials. Instead, at 100 K 1/Ti has a value 
that is ~ 10 3 times larger in A-(BETS) 2 FeCl 4 than in A- 
(BETS) 2 GaCl 4 . Futhermore, in A-(BETS) 2 FeCl 4 , l/T x 
has the T— dependence expected from the saturation of 
the magnetization of the Fe 3+ ions. 19 

These properties provide strong support to our inter- 
pretation that the dominant internal magnetic field at 
the protons in A-(BETS)2FeCl4 is the dipole field of the 
Fe 3+ ions, and that there is no anomalous splitting of the 
spectrum near 70 K. 

Other interpretations of a ferroelectric phase transi- 
tion at 70 K in A-(BETS) 2 FeCl 4 are based upon the 
discontinuity reported for its specific heat^ the division 
of the electron spin resonance g— factor from one into 
two values^ and the increase in the microwave dielec- 
tric constant £ It is not yet clear why the effect of such 
a transition on the magnetization of the Fe 3+ ions is so 
small that it is not seen in our 1 H-NMR spectrum mea- 
surements. 



V. CONCLUSIONS 

1 H-NMR spectrum measurements on a small (~ 4 /zg) 
single crystal of the organic conductor A-(BETS) 2 FeCl4 
using an applied magnetic field B = 8.9885 T parallel 
to the a-axis in the ac-plane over the temperature range 
2.0—180 K are reported. This work, and a preliminary 
report of itjiS are the first NMR reports that use a single 
crystal of A-(BETS) 2 FeCl 4 . The results provide the dis- 
tribution of the static local magnetic field at the proton 
sites in both the PM and the AFI phases. 

The experimental spectra have six main peaks and be- 
come progressively broadened and shifted as T is de- 



creased from 180 K to ~ 5 K. For T < T M i = 3.5 K 
(below the PM— AFI transition), an extra peak appears 
on the high frequency side, the details of the spectrum 
become smeared, and changes in the frequency shift and 
the rms linewidth are discontinuous, indicating a signifi- 
cant change in the static local magnetic field distribution 
at the proton sites on traversing the PM to AFI phase 
transition. 

The origin of the spectral features is attributed to the 
large dipolar field from the 3d Fe 3+ electron moments 
(spin Sd — 5/2, g w 2) at the proton sites. The main 
features of the spectra are successfully modeled with a 
mean field corrected Brillouin function. 

The value for J^d between each of the two nearest 
neighbor Fe 3+ ions obtained from this fit to the mea- 
surements is ~ — 0.85 K, which is close to the theoretical 
prediction*24 

No NMR evidence for an anomaly at 70 K reported 
earlier 7 on an aggregate of crystals is observed. 

It is suggested that the smaller features of the spectra 
that are not covered by this model are caused by the 
electron-electron interactions that are beyond the scope 
of this paper. 
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